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P { . ABSTRACT 

I ' We present a Spitzer Space Telescope spectroscopic survey of mass-losing carbon stars (and 

one oxygen-rich star) in the Large Magellanic Cloud. The stars represent the superwind phase 
> ■ on the Asymptotic Giant Branch, which forms a major source of dust for the interstellar 

medium in galaxies. The spectra cover the wavelength range 5-38 /im. They show varying 
combinations of dust continuum, dust emission features (SiC, MgS) and molecular absorp- 
tion bands (C2H2, HCN). A set of four narrow bands, dubbed the Manchester system, is used 
to define the infrared continuum for dusty carbon stars. The relations between the continuum 

• colours and the strength of the dust and molecular features are studied, and are compared to 

Galactic stars of similar colours. The circumstellar 7-{im C2H2 band is found to be stronger 
at lower metallicity, from a comparison of stars in the Galaxy, the LMC and the SMC. This 
is explained by dredge-up of carbon, causing higher C/O ratios at low metallicity (less O). A 
possible 10-^m absorption feature seen in our spectra may be due to C3. This band has also 
been identified with interstellar silicate or silicon-nitrite dust. We investigate the strength and 
central wavelength of the SiC and MgS dust bands as function of colour and metallicity. The 
line-to-continuum ratio of these bands shows some indication of being lower at low metallic- 
ity. The MgS band is only seen at dust temperatures below 600 K. We discuss the selection of 
carbon versus oxygen-rich AGB stars using the J— K vs. K— A colours, and show that these 
colours are relatively insensitive to chemical type. Metal-poor carbon stars form amorphous 
carbon dust from self-produced carbon. This type of dust forms more readily in the presence 
of a higher C/O ratio. Low metallicity carbon dust may contain a smaller fraction of SiC and 
MgS constituents, which do depend on metallicity. The formation efficiency of oxygen-rich 
dust depends more strongly on metallicity. We suggest that in lower-metallicity environments, 
the dust input into the Interstellar Medium by AGB stars is efficient but may be strongly bi- 
assed towards carbonaceous dust, as compared to the Galaxy. 

Key words: Galaxies: Magellanic Clouds; stars: stars: AGB and post- AGB - stars: carbon - 
stars: mass loss - stars: infrared dust 



1 INTRODUCTION 
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Low and intermediate-mass stars (M ~ 1-8 M©, hereinafter 
LIMS) make up >90% of all the stars which have died in the Uni- 
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verse up to the present time. During their late evolution, after enter- 
ing the high luminosity asymptotic giant branch (AGB) phase, the 
stars eject their hydrogen-rich outer layers during a phase of catas- 
trophic mass loss: the so-called superwind phase. After this final 
burst of activity, the star remains as a hot, compact white dwarf of 
mass 0.6-1.4 Mg. The expanding ejecta surrounding the star be- 
come ionized and form a planetary nebula, before dispersing into 
the interstellar medium. 

During the superwind, the star tends to be self-obscured by 
dust which forms in the ejecta. Dredge-up of primary carbon (pro- 
duced by triple-a burning) may turn the star into a carbon star. At 
the low effective temperature of AGB stars, the ejecta are largely 
molecular. The first and most stable molecule to form is CO: this 
locks up the least abundant of the C and O atoms. Other molecules 
and dust form from the remaining atoms. The chemistry changes 
dramatically as the C/ O ratio passes un i ty, and this affects both 
the gas-phase species jMillar et alj|200d: IWillacv & Millar|[T997l) 
and the dust composition (Treffers & Cohen 1974). In carbon stars, 
the C/O ratio (by number) is more than unity, and carbon-rich 
molecules (e.g. C2H2, HCN) and carbonaceous dust (Honv l2002h 
result, whilst in oxygen stars (C/0< 1), metal oxides and silicate 
dust <Camil20 021 form. 

The mass-loss process is important for three reasons. First, 
it determines the mass distribution of stellar remnants, includ- 
ing e.g. the lo wer mass limit of type-II supernovae progenitors 
(Ziilstra f2004h . Second, the physics and interplay of stellar pul- 
sation, shock waves, dust formation and radiation pressure which 
drives the mass loss are still little understood. Third, stellar mass 
loss drives galactic evolution through replenishment and enrich- 
ment of the Interstellar Medium (ISM). On this last point, mass- 
loss from LIMS contributes roughly half the total gas recycled by 
all stars (Maeder 1992), creates an amount of carbon roughly equal 
to that produced by super novae and Wolf-Rayet stars ( Dr ay et alJ 
l2003t iGavilan et alF2 005 ) and is the main source of carbonaceous 
interstellar dust ( Edm undsll200 it IPwetJI 1 998h . LEVIS are the only 
confirmed source of the primary nitrogen required to explain ob- 
served abundances in object s in th e early universe, are the main 
source of 13 C jHaiduk et alj|2005h and heavy s-process elements 
iGorielv & Mowlavi l2000|jWanajo_et_al. 2006), and are the major 
stellar source of lithium (Romano et al. 2001). 

The AGB and post- AGB evolution is dominated by the super- 
wind mass loss. Theoretical models (e.g. Sandin & H6fner H2003li 
cannot yet predict mass-loss rates from stellar parameters: instead 
observational knowledge is required of the mass-loss rates as a 
function of mass, luminosity and metallicity. However, the obser- 
vational relations are themselves not well calibrated: local Galac- 
tic stars have poorly known distances and unknown progenitor 
masses, while extra-galactic stars are in general too faint to de- 
tect the mas s-loss trace rs. Important work has been done on Bulge 
stars lOrtiz et aljEool . The main tracer for the superwind is the 
dust, but dust emits in the thermal infrared. Surveys by the Infrared 
Astronomical Satellite (IRAS) uncovered a number of luminous 
self-obscured AGB stars in the Magellanic Clouds lReid ll99ft 
Ziilst ra et alll996l) . A larger number of these stars was discovered 
more recent ly with th e Midcourse Space Experiment (MSX) satel- 
lite lEgan et alj200lh . The Infrared Space Observatory (ISO) made 
it possible to measure mass-loss rates, but only for the most mass ive 
stars and the highest mass-loss rat es (e.g.|va n Loon et al. 1999a). 

The Spitzer Space Telescope (Werner et al. 2004) for the first 
time provides the sensitivity necessary to obtain mass-loss rates for 
the whole range of AGB masses and luminosities in the Magellanic 
Clouds. These stars have well-known distances, meaning that abse- 
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Figure 2. The Mk, J — K diagram for the sources in th is paper. J an d 
K values are taken from Tables [Tlandl2l and from lvan Loon et al] 12005a ). 
Objects with lower limits only for the J— K colour are not included. Symbols 
indicate a source classification based on data described in this paper: open 
symbols: stars without a MgS feature, filled triangles: a weak SiC features 
and MgS, filled squares: strong SiC and MgS 

lute luminosities can be easily obtained. The general abundance of 
field stars may be estimated from the well-studied age-metallicity 
relations, and abundances have been derived explicitly for many 
star clusters. We have thus undertaken a survey of stars all along the 
AGB branch in the Large Magellanic Cloud and the Small Magel- 
lanic Cloud (hereinafter LMC and SMC respectively) with Spitzer. 
In this paper we present the LMC stars observed in our survey. Data 
for the SMC stars and the analysis of mass-loss rates will appear in 
subsequent papers. 



2 TARGET SELECTION 

The (A/k, J — K) diagram in Fig.0a) shows most of the AGB stars 
which are known mid-IR point sources in the LMC. Also shown 
is a representative sample of the much more numerous field stars 
without thick dust shells. The most numerous stars in any stellar 
population are low-mass stars. In Fig. 0a), these stars evolve up 
the AGB in M K at J - K < 2.5 and M K > -9, before they 
develop high mass-loss rates and thick dust shells and evolve to 
J — K ~ 7 and Mk ~ —5. Fig. 0a) also includes massive AGB 
stars (Al > 3 Mq , — 10.5 < Mk < —9), deriving from intermediate- 
mass progenitor stars. 

Stars observed in the current survey are shown in Fig. 0b). 
We selected targets in both the SMC and t he LMC from the Mid- 
course Space Experiment (MSX) catalog lEgan et aT]|200lll . the 
IRAS point source catalog, the ISO catalog of variable AGB stars 
in the SMC ICioni e Ud]|2Q03l), long-period variables in the Mag- 
ellanic Clouds (e.g. IWood et al"]ll983t IPeast et aljfl989l) . and the 
intermediate-age clusters NGC419 and NGC1978, both of which 
have a known mid-infrared source. Targets were selected to span 
the long sequence of increasing mass loss in LIMS which stretches 
from (J-K,M K ) ~ (2,-9) to (7,-4) in Fig. 0a). (The bright 
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Figure 1. (a): The Mr, J-K diagram for a sample of known mid-IR sources (MSX and IRAS) in the LMC (large points), together with field stars (small 
symbols) from an area ~ 60 times smaller than the area from which the mid-IR sources were selected. AGB stars are confined approximately to the region 
below the line marked 'AGB limit". Stars above this limit are foreground stars and supergiants: the supergiant stars have masses M > 8Mq. Spectral types 
were those available before the Spitzer observations were made (b): similar to (a), except that te objects shown are the point sources actually observed, and 
the spectral types come from recent ground-based and Spitzer spectra. In both panels, distance modulii of 18.5 and 18.9 have been assumed for the LMC and 
SMC, respectively. J and K photometry is from various sources existing prior to the Spitzer observations. 



O-rich object near (J-K,M K ) ~ (1,-10.7) in Fig. [3b) was 
not a target but is an object which was accidently observed when 
Spitzer selected it in the peak-up field). It turns out that only two 
of the LMC sources were oxygen rich, one of wh ich was found to 
be a massive young stellar object ( van Loonetall2 Q05b). Here we 
present the carbon-rich sources in the LMC. The SMC sample will 
be presented elsewhere. 

TableQlists the 29 stars included here. The objects have been 

included in various catalogues over time, accumulating different 

names. We adopt the oldest available name for each object, but Ta- 
rn I . 1 1 I I I 
bleHJgives the cross references to names used by Reid et al. 1 1990), 

lEean et all 1200 ll) and in the MSX catalogue. The adopted positions 
are from 2MASS counterparts and differ slightly from the MSX po- 
sitions. 

The JHKs photometry in Table Q is taken from the 2MASS 
All-sky catalog. It is based on a single epoch only and does not 
sample the variability. Where further data is available, we list JHKL 
photometry averaged over t he cycle in Table[2] Als o listed are pul - 
sation periods taken from Whiteloc k et alj j2003l).IWood | l998|) . 
iNishida et all feOOOl) . Ilta et alj 120041) and iGroenewegerj |20o4) . 
Where available, periods from Whitel ock et al] 12003) are listed in 
preference because of the larger number of observations used. Pho- 



tometry is taken from the first three sources, with values from Wood 
listed in preference, where available. Especially the L-band filters 
differ substantially between the sources; the photometry is in the 
original system. For MSX LMC 967, the photometry is an average 
of two observa tions only: 2MASS and a single measurement from 
Ilta etaljExil . 

Fig. |U shows the Mk versus J — K colour-magnitude dia- 
gram, using pho tometry listed above supplemented with data from 
Ivan Loon et alJ l2006), The stars follow the line of increasing mass 
loss defined in Fig.Qa). Neither figure includes the stars which are 
faintest at K, as these show only uper limits for J. 



2.1 Comments on individual stars 

Three stars in TableQare located in the cluster NGC 1978, as in- 
dicated by their names. The cluster was studied by iTanabe et alj 
( 1997); it has an age of 2.5 10 9 yr, a turn-o ff mass of about 1.5 Mr 
and a metallicity [Fe/H] = -0.66 ± 0.22 Jvan Loon et aJ2005al) . 
One star, MIR-1, is a known obscured LPV , one other is an MSX- 
detec ted star, and the third (IR4, van Loon et al . 2005al lWilletai] 
1995) is an AGB star with no listed MSX point-sou rce counterpart. 
This last star is seen on MSX images, from which Ivan Loon et alj 
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Table 1. Observed LMC targets: names, adopted coordinates, and photometry. The MSX photometry is for band A (8.2 micron) and refers to version 2.3 of 
the catalogue. The zeropoint for t he MSX band A is ta ken as 58.46 Jv {Egan et all200ll) . JHK is taken from 2MASS, except for MSX LMC 494 and NGC1978 
MIR1 where the values are from llta et all2004l) and Ivan Loon et all2005al) respectively. M bo i (Section [3.H assumes a distance modulus of 18.5. 



Adopted name 


MSX source 


other name 


RA Dec 
(J2000) 


J 

mag 


H 

mag 


Ks 
mag 


A 
mJy 


A^boi 


carbon rich 






















MSX LMC 219 


051119.5- 


-684227 




05 11 19.52 


-68 42 27.9 


>15.61 


> 14.67 


13.14 


204 


-4.76 


IRAS 05132-6941 


051250.8- 


-693749 


MSX LMC 223 


05 12 51.07 


-69 37 50.3 


>17.81 


>15.68 


13.36 


259 


-5.01 


MSX LMC 349 


051726.9- 


-685458 




05 17 26.94 


-68 54 58.2 


> 17.44 


> 16.88 


14.82 


193 


-4.97 


MSX LMC 341 


052100.5- 


-692054 




05 21 00.37 


-69 20 55.3 


> 17.49 


15.9: 


13.151 


210 


-4.74 


MSX LMC 663 


052244.5- 


-693826 




05 22 44.00 


-69 38 28.1 


12.60 


11.33 


10.20 


75 


-5.40 


MSX LMC 494 


052309.6- 


-691744 




05 23 09.11 


-69 17 49.1 


13.08 


11.71 


10.51 


94 


-4.98 


MSX LMC 441 


052438.7- 


-702357 




05 24 38.62 


-70 23 57.1 


>17.16 


> 16.44 


14.48 


200 


-4.92 


MSX LMC 443 


052505.9- 


-701011 




05 25 05.69 


-70 10 10.6 


> 14.54 


13.72 


11.54 


150 


-4.68 


MSX LMC 601 


052650.9- 


-693136 




05 26 50.83 


-69 31 36.9 


13.88 


11.99 


10.52 


172 


-5.14 


NGC 1978 IR4 






WBT 815 


05 28 44.50 


-66 14 04.0 


11.49: 


10.34: 


9.68 


70: 


-5.35 


NGC 1978 MIR1 
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16.67 




12.45 


124 


-4.75 


NGC 1978 IR1 






WBT 1268 


05 28 40.17 


-66 13 54.2 


>13.78 


>12.78 


11.73 


77 


-5.03 


IRAS 05278-6942 


052724.3- 


-693944 


MSX LMC 635 


05 27 24.12 


—69 39 45.0 


>17.53 


>15.59 


12.35 


683 


-6.57 


MSX LMC 754 


052811.4- 


-703359 




05 28 11.48 


-70 33 58.7 


15.98: 


13.70 


11.89 


93: 


-4.91 


MSX LMC 679 


052848.5- 


-694801 




05 28 48.62 


-69 48 01.3 


>15.84 


14.57 


12.48 


121 


-4.93 


MSX LMC 743 


053454.1- 


-702925 




05 34 53.74 


-70 29 24.8 


>14.13 


> 13.22 


13.27 


165 


-5.27 


MSX LMC 749 


053527.1- 


-695229 




05 35 26.86 


-69 52 27.9 


15.06 


12.93 


11.05 


142 


-5.56 


MSX LMC 967 


053637.2- 


-694725 




05 36 36.71 


-69 47 22.6 


13.48 


12.00 


10.67 


87: 


-5.45 


IRAS 04557-6753 






MSX LMC 1238 


04 55 38.98 


-67 49 10.7 


>16.18 


14.49 


12.397 


157 


-5.05 


IRAS 05009-6616 






MSX LMC 1278 


05 01 04.43 


-66 12 40.4 


> 15.64 


14.57 


12.40 


240 


-5.57 


IRAS 05291-6700 






GRRV 38 


05 29 07.60 


-66 58 15.0 


12.45 


10.91 


9.91 


81: 


-4.91 


IRAS 05295-7121 






MSX LMC 692 


05 28 46.62 


-71 19 12.5 


16.80: 


14.32 


12.18 


204 


-5.24 


TRM88 






MSX LMC 310 


05 20 19.38 


-66 35 47.8 


14.76 


12.81 


11.10 


169 


-5.15 


IRAS 05112-6755 






MSX LMC 44, TRM4 


05 1 1 10.47 


-67 52 10.5 


16.41: 


14.07 


11.69 


451 


-5.59 


IRAS 05190-6748 






MSX LMC 307, TRM20 


05 18 56.26 


-67 45 04.4 


>18.21 


>15.91 


13.10 


209 


-5.61 


IRAS 05113-6739 






MSX LMC 47, TRM24 


05 11 13.89 


-67 36 16.1 


>17.72 


14.77 


12.49 


233 


-5.18 


TRM72 






MSX LMC 29 


05 11 38.65 


-66 51 09.8 


14.34 


12.36 


10.74 


122 


-5.43 


IRAS 05360-6648 






MSX LMC 872, TRM77 


05 36 01.24 


-66 46 39.7 


> 18.04 


16.07: 


13.28 


192 


-4.73 


oxygen rich 

IRAS 05003-6712 






MSX LMC 1280 


05 00 19.00 


-67 07 58.0 


12.04 


10.46 


9.32 


205 


-5.42 



<2005d) estimate a flux density of 70 mJy in MSX band A (6.8-10.8 
micron). 

IRAS 04557-6753 is also located in a ( small) cluster, KMHK 
285 <Bica et all 19991: Ivan Loon et all2 005a). but this cluster is lit- 
tle studied. The age is estimated as > 10 9 yr and the metallicity is 
not known. 

MSX LMC 494 is identified with a large amplitude vari- 
able observed with the Optical Gravitational Lensing Experiment, 
OGLE 052309.18-691747.0. 

LMC MSX 663 is also a l ong-period OGLE v ariable but its 
OGLE amplitude is very small iGroeneweg en 2004). I t is an opti- 
cally bright AGB star previously classified as an S star ICioni et all 
2001). However, optical spectra show that it is a symbiotic carbon 
star with strong emission lines (Wood et al, in preparation). 

IRAS 0529 1—6700 was id entified with a nearby variable star, 
GRRV 38 iGlass & Rei3ll985l) . but the coordinates of the IRAS 
source differ by 30 arcsec from this star. The IRAS point source is 
likely unrelated. There is confusion with an extended mid-infrared 
source approximately 6 arcminutes east. MSX only has the stel- 
lar source, but nothing at the position of the IRAS source. The 
identification of the star I Zijlstra et al. 1996) most likely occurred 
serendipitously. GRRV 38 is a carbon star Ivan Loon et alll999bl) : 



the infrared colours indicate little or no obscuration. Below we will 
use the name GR RV 38 for this obje ct. 

For TRM 88. IWoodil99i) and Whitelo ckeTaT] 120031) report 
main periods of 565 and 544 days res pectively. T he light curve is 
affected by long-term secular chan ges j Whitelo ck et all2003l) due 
to variable circumstellar extinction iFeast et^lT2003l) . 

IRAS 05009-6616 and TRM 72 also show evidence for long- 
term variations in the light curves jWhitelock et all2003l) . 

TRM 72 should not be confused with the well-studied CH star 
WORC 106 iFeast & Whitelocklll992l) which is located 30 arcsec 
away. 



3 OBSERVATIONS 

We observed our targets with the Infrared Spectrograph (IRS) 
iHouck et all2004l) on the Spitzer Space Telescope using the low- 
resolution modules Short-Low (SL) and Long-Low (LL), which 
together cover wavelengths from 5 to 38 /im. The slits in SL are 
roughly perpendicular to LL, and each module has two separate 
apertures, with the longer wavelengths observed in first order and 
the lower in second order. In order of increasing wavelength, the 
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Table 2. Ph ase-averaged photometry and periods. Periods are taken 
from IWhitelock et all 120031) . IWood Il998l) . iNishida et all l2000t) . and 
iGroenewegenl <2004l) . 



Adopted name 


J 


H 


K 


L 


P 




mag 


mag 


mag 


mag 


days 


IVIOYY I J L V 1 V . UUJ 










455" 


MSX LMC 494 










458 


MSXLMC 601 










546 


NGC 1978 IR1 


14.26 


12.01 


10.51 




491 


MSX LMC 749 










329, 589 


MSX LMC 967 


(14.2) 


(12.5) 


(10.8) 




572 


IRAS 04557-6753 




13.97 


11.78 


9.11 


765 


IRAS 05009-6616 




13.06 


11.29 


8.94 


658 


GRRV 38 


13.04 


11.33 


10.17 


8.86 


483 


IRAS 05295-7121 




12.96 


10.99 


8.83 


682 


TRM88 


14.22 


12.10 


10.30 


8.26 


544 


IRAS 05112-6755 


17.65 


15.11 


12.29 


8.89 


830 


IRAS 05190-6748 




16.11 


12.72 


8.71 


939 


IRAS 05113-6739 


16.84 


14.63 


12.17 


8.93 


700 


TRM72 


15.50 


13.33 


11.39 


8.92 


631 


IRAS 05360-6648 




15.50 


12.87 


9.24 


538 


IRAS 05003-6712 


12.90 


11.26 


9.95 


8.49 


883 



0/2 



^ 0.1 



four spectral segments are known as SL2, SL1, LL2, and LL1. In 
addition, a "bonus" order covers the overlap between the first and 
second orders in each module; this order is a small piece of the first- 
order spectrum obtained when the source is in the second-order 
aperture. When the source is in one aperture, the other aperture is 
exposed to a section of nearby sky, providing a background mea- 
surement, e.g. in SL2 when the source is in the SL1 aperture. Each 
source was observed in two nod positions in each aperture. The 
resolution varies between 64 and 128, depending on the order and 
wavelength. The slit dimensions are approximately 3.5" x 57" for 
SL and 10.5" x 168" forLL (width x length). 

Different exposure times were used depending on expected 
target flux. For SL, 12 to 120-sec exposure times were used per 
segment. For LL1, exposure times ranged from 60 to 960 sec, and 
for LL2 from 180 to 960 sec. Observations were taken in standard 
staring mode. 

Flat-fielded images were generated by the Spitzer Science 
Center (SSC). We replaced bad pixels with values estimated from 
neighbouring pixels (using the imclean.pro package), and re- 
moved the background emission by differencing images aperture- 
by-aperture in SL (e.g. SL1— SL2 and vice versa), and nod-by-nod 
in LL. Spectra were extracted from the images using the software 
provided with the Spitzer IRS Custom Extractor (SPICE) available 
from the SSC. They were calibrated making a full spectral cor- 
rection using HR 6348 (K0 III) as a standard star in SL and HR 
6348, HD 166780 (K4 III), and HD 173511 (K5 III) in LL. Indi- 
vidual images were extracted separately and then co-added. Nods 
were combined, then spectral segments were combined by making 
scalar multiplicative corrections for pointing-induced throughput 
errors (normalizing to the best centered segment). The accuracy of 
the overall spectral shape is limited to a few percent. The match- 
ing across the different sections (e.g. LL2 to LL1) may leave some 
residuals. 

The standard wavelength calibration is accurate to 0.06 micron 
in SL and 0. 15 micron in LL. Flux calibration is accurate to ^ 5 per 
cent. The response curve is known to a few per cent, but residuals 
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Figure 5. Spectrum of the symbiotic C star MSX LMC 663 



at this level may remain. The largest residual calibration errors are 
due to pointing uncertainties. The LL array is affected by fringing 
which in most cases has been corrected. More details on calibration 
issues can be found in lSloan et alJ 1200 3l) . 

The spectra of the carbon stars are shown in Fig.|3|and |4] They 
are ordered by dust temperature, from blue spectra to red. The spec- 
tra show numerous absorption and emission bands of molecules 
and dust resonances, superposed on the dust continuum. We will 
first discuss these bands, separated into dust and molecular fea- 
tures. The following section will discuss the continuum spectral 
energy distribution. 

Some objects show residual interstellar lines, due to in- 
complete cancellation of the background emission. These lines 
are found at 33.5 /im ([S III]), 34.8 fim ([Sill]), and 36.0 /im 
([Nelll]). An apparent feature at 30.5 /im seen in a few spectra 
is an artifact, due to a hot pixel. 

The object MSX LMC 663 exhibits an IRS spectrum which is 
very different from the others, as shown in Fig. [5] It is affected by 
interstellar emission in the background measurement, causing the 
offset measurements to disagree; however, the rise beyond 25 /J,m 
appears to be intrinsic to the star and not due to background emis- 
sion. Shortward, the spectrum shows little evidence for dust but the 
absorption bands are different in shape from those of the other car- 
bon stars. 

IRAS 05003-6712 is the only O-rich AGB star in our sample 
(one other object was classified as a massiv e young stellar object, 
and is published in Ivan Loonetai][2005b). It has an OH maser 
iMarshall et ai] f2004 ) but only the blue-shifted component is de- 
tected. The Spitzer spectrum is dominated by the silicate emission 
bands at 10 and 18 /im, with perhaps a hint of CO2 absorption at 
13.5 fim and SiO absorption at 8.6 /im . (SiO bands tend to be weak 
in LMC Miras: iMatsuura et alj2005T) . The l0-p,m feature is typical 
for a lower mass-loss rat e: it is com parable to that of the LMC AGB 
star IRAS 04544-6849 iDiikstra et all2005h with a mass-loss rate 
of M = 7 x 10" 6 M vr" 1 Ivan Loon et all 19993) . . 



3.1 Bolometric magnitudes 

For a number of our targets, bolometric magnitudes were derived 
by IWood ll998t) and IWhitelock etal] <2003l) . Table □ lists the 
bolometric magnitudes derived from the Spitzer spectra and JHKs 
2MASS photometry. Where multi-epoch photometry was available 
(Table|2j, those were used instead of the single-epoch 2MASS data. 
The spectra were smoothed by a factor of 20. The J-band upper 
limit was used as a real magnitude where necessary: this has no 
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3. Spectra of the 15 bluer carbon-rich AGB stars observed in the LMC. The top five lines show the model molecular spectra from ljorgensen etHl] 
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Figure 4. Spectra of 1 1 redder carbon-rich AGB stars observed in the LMC, and spectra of the blue ('naked') carbon star NGC 1978 IR4 which was observed 
twice. The top five lines show the model molecular spectra from lj0rg ensen etlifl bood) 
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Figure 6. Spectrum of the oxygen-rich AGB star IRAS 05003-6712 
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Figure 7. Distribution of bolometric magnitudes 



4 EMISSION AND ABSORPTION BANDS 
4.1 Dust bands 

Carbon stars show two major dust bands in this spectral region, 
at 11.3 and 30 /im. Both features are seen in the majority of the 
current sample, although not in all. Two other spectral features have 
been suggested as being due to dust, but below we argue against this 
possibility. 



4.1.1 SiC dust emission 

The w ell-known 1 1.3-/xm band is du e to solid a-SiC (e.g. 
IPapoular et alll998tlcTement et all2 003). The central wavelength 
shows some variability in Galactic carbon stars, ranging between 
11.15 and 11.7 fan. Longer wavelengt hs are correl a ted with lower 
dust (c ontinuum) temperatures (Bar on et alJll987j) . Iciement et alJ 
(2003) find that the wavelength shift to the red can be reproduced 
in SiC particle conglomerates, and they suggest that conglomerates 
form at higher density in the wind. In this case the correlation with 
the dust temperature is an indirect one, caused by a larger optical 
de pth in the dust she ll. SiC absorption components are discussed 
bv lSpeck etaLh2005h . 



4. 1.2 MgS dust emission 

The 30-/itn resonance is commonly seen in carbon-rich envi- 
ronments, in AGB stars, post-AG B stars and planetary nebulae 
iForrest et alJ 1 1981 1 : [Honvl |2002|> . It consists of two subpeaks, 
at 26 and 33 /im (Volketal. 2002). The two peaks always ap- 
pear together, arguing for an origin in the same miner al. The 
resonance is attributed to MgS (Goebel & Mos ele^ll985[) which 
forms via grain-surface reaction between Mg and H?S~ lNufh et alJ 
ll985j); the two pe aks have been attributed to different grain shapes 
iHonv et all2002l) . A contribution from CaS to the features cannot 
be fully excluded, but the absenc e of any feature in oxygen-rich 
environments argues against this I Goebel & Moselev 1985J). The 
peak wavelength of the MgS feature shifts with temperature, which 
make s it an important tem perature diagnostic iHonv & Bouwmanl 
120041) . IHonv et all |2002) argue that solid MgS may survive into 
the ISM. 



measurable effect on the result (changing it by 1 mag affected the 
bolometric flux by less than 1 per cent). Zero flux was assumed at 
frequencies of Hz and 3 10 16 Hz. A fourth-order polynomial was 
fitted (using IDL) to the spectra and data points and the total flux 
was obtained by integration under the fit. For the blue carbon stars, 
an uncertainty is introduced by the lack of data points between 2 
and 5/j.m. Variability is the main limitation to the accuracy: the 
K-ban d magnitude can vary by up to 1-2 mag in ex treme cases 
<Le Bertrdll99i IWood[l99i Iwhitelock et alll2003l) . The result 
was converted to absolute magnitudes assuming a distance mod- 
ulus of 18.5 (e.g., Alves et al. 2002). The distance modulus varies 
by ap proximately 0. 1 magnitude over the face of the bar lLah et all 
2005) but we did not correct for such geometric effects. 

The histogram of the bolometric magnitudes is shown in Fig. 

□ 



4.1.3 Doubtful features 

Two objects show indications of a broad, weak hump between 16 
and 23 ^m, with unclear identification. It is most noticeable in the 
symbiotic C star LMC MSX 663, and may be seen weakly in MSX 
LMC 754, a carbon star with hot dust. Neither object shows evi- 
dence for MgS. The feature is likely an artifact caused by the broad 
molecular C2H2 absorption band to the blue. 

A possible emission feature at 8.6 /im has been described in 
IRAS LRS (Low Resolution Spectrometer) spectra of Galactic car- 
bon stars iWillemsil98"8llSloan et alll998h . The feature correlates 
in strength with the SiC band and the 3 fim b and. However, its re- 
ality is disputed: strong absorption at 7.5 /im lAoki et alll99 9) to- 
gether with 10 /^m absorption can mimic an emission band. The 
spectra in our LMC sample show strong absorption at 7.5 fim and 
absorption at 10 fim, making it more likely that the 8.6 fim peak is 
simply the continuum between two molecular bands, which we will 
assume below. 



4.1.4 10 micron absorption 

This absorption feature, which has been attributed to a dust compo- 
nent, is discussed in Section l4~3l 



4.2 Molecular bands 

The molecul ar bands are all at the short wavelength range of the 
IRS spectra. iGautschv-Loidl et alJ 120041) discuss the contribution 
of different molecules to the 5-20-^im wavelengt h range. Exam- 
ple mo dels for individual bands were calculated by J0rgensen et al. 
(2000): their models are shown at the of Figs.[3]and|4] 



4.2.1 13.7 micron 

The narrow 13.7-/^m band is prominent in most of the s pectra. 
This band i s due to th e Q-branch C2H2 V5 transition iTsuii 1984; 
lAoki et alJll999l ICernicharo etafl ll999). The narrow band is sit- 
ting in the much broader P,R branches, giving rise to an absorp- 
tion feature extending from 12 to 16-/nm, which is best seen in 
the spectrum of IRAS 05190-6748. The 11.3-/xm SiC emission 
feature makes it difficult to judge how stro ng this broad a bsorp- 
tion b and is. Models of Galactic carbon stars iGautschv-Loidl et alJ 
120041) predict a much stronger broad band than is observed. This 
can be explained if the mass-losing envelope adds a molecular 
emission feature. The sharp 13.7-/im band is very strong in our 
spec tra, while in Galactic stars i t has been dubbed the 'weak' fea- 
ture iGautschv-Loidl et alj20 04). The sharp feature originates from 
colder gas, while hot (up to photospheric temperatures) gas mainly 
shows the broad band lJ0rgensen et al. 2000). 



4.2.2 14.3 micron 

The 14.3-/xm HCN band seen in Galactic carbon stars lAoki et ail 
Il999t ICernicharo et alj[l999l) is weak in our sample. It is proba- 
bly present in the spectrum of TRM 72. Otherwise the feature is 
not clear, and is certainly much weaker than the C2H2 band. This 
is similar to the situation in Galactic stars, w here its weakness is 
attributed to in-fill by an emission component lAoki et alll 999). 



4.2.3 7.5 micron 

A noticeable absorption band is present at 7.5 fim. Galactic car- 
bon stars have contributions in this wavelength region from C2H2, 
H CN and CS iGoebel et alJll98lL lAoki et aljll99l) . The sample 
of lAoki et alJ i^9^T consists of blue carbon stars, which differ in 
effective temperature from our sample. 

The shape of the two shallow features found in our spectra 
resembles C2H2. The diatomic molecule CS shows a band head 
which we do not see, and the HCN band is located at shorter 
wavelength, typically centered around 7.0 (although this central 
wavelength depends on the excitation temperature). Comparison of 
our LMC spectra and Galactic spectra indicates a difference in the 
shape of the 7.5-/im band, supporting different origin s for Galactic 
and LMC carbon stars. iMatsuura et alJ i2002alEo05T) show that at 
low metallicity, C2H2 becomes more abundant compared to HCN. 
Following this, we identify the observed band for LMC stars with 
C2H2. 
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Table 3. Wavelengths used to estimate the continua for the SiC and C2H2 
spectral features. 



Features 


A [/jm] 


Blue continuum [/tm] 


Red continuum [(im] 


C 2 H 2 


7.5 


6.08-6.77 


8.22-8.55 


SiC 


11.3 


9.50-10.10 


12.80-13.40 


C2H2 


13.7 


12.80-13.40 


14.10-14.70 



4.2.4 5 micron 

At the blue edge of the spectra the spectral energy distributions 
drop sharply, due to a strong absorption band which extends beyond 
the blue limit of the wavelength range. There are two molecules 
which contri bute at this wa velength in carbon-rich environments: 
CO and C3 1 J0rgensen et al. 2000). CO is of photospheric origin 
and is present even in K stars. Its presence in e.g. the unobscured 
source GRRV 38 does not require a circumstellar shell. C3 also is 
photospheric in carbon stars. The two bands differ at the red end: 
C3 ends at 6 /im, while CO tapers off slowly up to 8 /im. 

SiS has three bands at 6.6-6.7 fim but the IRS resolution is not 
high enough to detect these; R = 2000 would be required. 



4.3 The 10-micron band 

An absorption feature at 10 /im may be present in most stars. 
It is s een o nly occasionally in Galactic carbon stars IVolk et alJ 
1200(1 ISloan et aljll993l . I J0rgensen et ail §000) find it in R Scl 
but no t in three other carb on stars they study. The origin is dis- 
puted. Clement et al. 1 2005) fit this band us ing solid silico n-nitrite 
particles. SkNd. lSpeck etaflll997l) and lVolk et aH i2000) suggest 
interstellar silicate absorpti on; to fit the shape, additional SiC ab- 
sorption is required. J0rge nsen et alJ J2OOO) identify the band in- 
stead with C3 (see their Fig. 11). 

The presence of the band in our LMC stars argues against the 
interstellar origin proposed by Sp eck et al]ll997l) . because the line- 
of-sight extinction to the LMC is too small to give significant sil- 
icate absorption. A contribution from SiC absorption, as required 
by the silicate hypothesis, is also ruled out: this would only be in 
absorption in dense, cool shells, while we see the 10-/im absorption 
even in blue, unobscured stars. 

This last point, the presence of the feature in all stars, argues 
against a dust source an d bring into question the suggested identifi- 
cation iQement et alj2005l) with silicon-nitrite. A molecular origin 
is favoured, in particular a molecule with a photospheric origin. 

We finally note the prevalence of the feature in our sample, as 
compared to Galactic stars. This points towards a molecule which is 
more abundant in the L MC (at low metallicity) than in the Galaxy. 
Matsuura et al. ( 2002a) have shown that this is the case for carbon- 
bearing molecules, due to a higher C/O ratio. The combined effect, 
at lower metallicity, of increased efficiency of third dredge-up of 
carbon 1 Wood ! 19811) and lower oxygen abundance leads to higher 
C/O ratios at lower metallicity. 

The identification suggested bv lj0rgensen et alH20 00) fits our 
available constraints. We therefore favour C3 as the probable car- 
rier. However, confirmation will require a laboratory measurement 
of its 10-/im spectrum. 
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Figure 8. The continuum underlying the MgS feature is subtracted using 
a black body, whose temperature is defined using the [16.5] — [21.5] color. 
The hatched bars display the wavelengths used to determine the [16.5] and 
[21.5] colors. The dashed line represents the black body fit to the spectrum 
of IRAS 04557-6753 taken as an example. 



Figure 9. The [6.4] — [9.3] versus [16.5] — [21.5] colour-colour diagram. 
Symbols are defined in Fig.|2| The drawn line represents a series of Dusty 
models, shifted down by 0.15 mag in [16.5]— [21.5]. 



Table 5. Mid-infrared continuum bands for carbon stars for the so-called 
'Manchester system'. The last column gives the adopted flux density corre- 
sponding to zero magnitude. 



central A [/im] 


A-range [fim] 


Fo [Jy] 


6.4 


6.25-6.55 


96.5 


9.3 


9.1-9.5 


45.7 


16.5 


16-17 


15.4 


21.5 


21-22 


9.1 



4.4 Band strengths 

We define the strength of the molecular absorption bands through 
their equivalent widths, and quantify the dust emission bands by 
their integrated line-to-continuum ratio. We use the adjacent con- 
tinuum to define a linear (sloping) continuum distribution across 
each band and define the equivalent width in the usual way. The 
wavelenghts used to estimate the continuum for the C2H2 (7.5 /Ltm), 
SiC (1 1.3 /im) and C 2 H 2 (13.7 (jxn) are listed in Tabled For the 
13.7-/im band, the equivalent width mainly refers to the narrow 
component: the broad component extends beyond the chosen con- 
tinuum wavelengths. 

The continuum underneath the MgS band cannot be measured 
in this way, because the band is too broad and extends beyond 
the red edge of the spectral coverage of the IRS. We therefore in- 
stead use a black body to approximate the continuum, and give the 
strength as a line-to-continuum ratio integrated over the feature up 
to 38 /im. The black body is defined using the [16.5] — [21.5] colour 
discussed below (see Fig.[8j. Table|4]lists the resulting values. 



5 CONTINUUM DEFINITION: THE MANCHESTER 
SYSTEM 

Much of the observed spectral range is covered by the molecu- 
lar and dust bands discussed above. This leaves the continuum 
poorly defined: one needs to select specific wavelength regions 
which avoid these bands. The broadband IRAS filters are clearly 
not a good approximation. The blue continuum is most affected by 
molecular bands; the red continuum is in many cases obliterated 
beyond 15 /J.m by the 30-/im MgS band. 

We have selected four narrow regions which can be used for 
dusty carbon stars (not for oxygen-rich stars!). The regions avoid 
most, but not all, of the obvious absorption and emission bands. The 
continuum bands are defined in Table [5] The bands will be called 
the Manchester system. Table|5|also lists the adopted zeropoint for 
the bands. These are extrapolated from the standard bands, and aim 
to return a zero colour for a Rayleigh- Jeans tail. 

The 6.4-/im band is situated between the SiS bands at 6.6 /im 
and CS at 7 um on the red side, and the CO band at 5 fim on the blue 
I Aoki et al. 1998). The wing of the CO feature may still affect the 
chosen band. There is also overlap with the PAH 6.2-fim band but 
this band requires UV excitation and is seen in planetary nebulae 
and post-AGB stars but not in AGB star spectra. 

The 9.3-/im band is blueward of the SiC dust band, and red of 
the CS bands and C 2 H 2 bands at 7-8 /jm lAoki et alll999l) . There 
are also bands from molecules such as NH3 and C 2 H 4 which ma y 
depress the continuum by a few per cent ICemicharo et all2001al) . 
However, it is situated in the C3 band discussed above, which is 
the main uncertainty in this spectral region and is significant for 
spectra dominated by photospheric emission. For these stars, the 
9.3-/im band will underestimate the true continuum. 

The two bands at 16.5 and 21.5 /im do not overlap with any 
strong molecular or dust features in carbon-rich material: it is sit- 
uated just redward of the broad P-branch 14-micron C 2 H 2 band, 
which extends to 16 /im. There is a weak hydrocarbon feature at 
the blue edge of the chosen band ICemicharo et aljboOlbt) . The 
strong 21-/im dust feature which is still unidentified is seen in post- 
AGB stars only. The Manchester system should therefore be used 
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Table 4. Strength of the molecular and dust features, in terms of either the equivalent width (EW) for molecular absorption bands, and the line-to-continuum 
ratio (L/C) for the dust emission bands. The central wavelength of the SiC band is also listed. The final column gives the continuum (black-body) temperature, 
derived from the [16.5] — [21.5] colour listed in Tablelol 



target 


EW (7.5 pm) 


EW(13.7/xm) 


L/C (SiC) 


A c 


L/C (MgS) 


T(K) 


TRM 72 


f)S2+0 005 


068+0 006 


n 071+0 002 

vj. yj i *j _i_v/. \j\JsL 


1 1.37 


o 


1018 


IRAS 04557 


0.156±0.004 


0.066+0.003 


0.123+0.003 


11.30 


0.412 


426 


IRAS 05009 


0.069+0.005 


0.062+0.004 


0.159+0.003 


1 1.31 


0.249 


619 


IRAS 05 1 1 2 


0.121+0.005 


0.066+0.003 


0.099+0.002 


11.34 


0.271 


396 


IRAS 05113 


0.133+0.008 


0.049+0.004 


0.094+0.004 


11.34 


0.215 


540 


IRAS 05132 


0.092+0.006 


0.071+0.005 


0.117+0.004 


11.27 


0.198 


557 


TR AS 05190 


n i 004 


074+0 004 


1 04+0 003 


1 1.37 


0.374 


398 


IRAS 05278 


0.082+0.004 


0.075+0.003 


0.111+0.003 


11.32 


0.407 


421 


GRRV 38 

V J 1V1V V JO 


n 91 8+0 023 


1 1 s+o 025 


019 

\J . Itj 1 \_J.\J 1 7 


11.14 


o 


858 


TR AS 05795 


n 0Q2+0 007 


057+0 004 


086+0 007 


1 1.34 


0.162 


462 


TR AS 05360 


n Q92+0 004 


05 1 +0 006 


1 80+0 006 


11.19 


0.624 


395 


MSX LMC 219 


0.111+0.005 


0.061+0.004 


0.108+0.003 


11.29 


0.336 


393 


MSXLMC341 


0.083+0.007 


0.063+0.005 


0.153+0.004 


11.16 


0.637 


482 


MSX LMC 349 


0.079+0.007 


0.049+0.005 


0.173+0.004 


11.18 


0.713 


470 


MSX LMC 441 


0.129+0.007 


0.056+0.004 


0.125+0.003 


11.26 


0.446 


427 


MSX LMC 443 


0.135+0.008 


0.054+0.007 


0.184+0.006 


11.18 


0.459 


467 


MSX LMC 494 


0.109+0.007 


0.039+0.003 


0.080+0.003 


11.31 





780 


MSX LMC 601 


0.016+0.005 


0.060+0.005 


0.056+0.003 


11.30 





703 


MSX LMC 679 


0.122+0.006 


0.044+0.004 


0.151+0.004 


11.22 


0.477 


415 


MSX LMC 743 


0.158+0.005 


0.042+0.004 


0.156+0.004 


11.24 


0.383 


455 


MSX LMC 749 


0.120+0.003 


0.076+0.003 


0.058+0.004 


11.25 





865 


MSX LMC 754 


0.121+0.006 


0.053+0.004 


0.146+0.006 


11.32 





726 


MSX LMC 967 


0.085+0.003 


0.058+0.003 


0.086+0.003 


11.30 





742 


NGC 1978 MIR1 


0.157+0.004 


0.064+0.006 


0.162+0.005 


11.30 


0.528 


446 


NGC 1978 IR1 


0.123+0.003 


0.070+0.004 


0.094+0.004 


11.23 





1179 


NGC 1978 IR4 


0.172+0.010 


0.082+0.009 


0.107+0.006 


11.50 





>10 3 


TRM 88 


0.051+0.006 


0.064+0.005 


0.065+0.004 


11.26 





684 



for carbon-rich AGB stars only. The regions also avoid potential 
instrumental problems in the IRS, such as the LL1/LL2 interface. 

The continuum fluxes of all observed targets in the Manch- 
ester system are listed in Table [6] These are derived by integrat- 
ing the IRS spectra over the bands. We give the colours using the 
ze ropoints of Table |5] Similar colours for SMC stars are given 
in lSloan et alj 12006). Table |6] lists a black-body colour tempera- 
ture derived from the [16.5] — [21.5] colour. For the bluest spec- 
tra, the values become poorly defined as the red colour approaches 
the Rayleigh-Jeans limit, but for the majority of objects a rea- 
sonable value is derived. Apart from NGC 1978 IR1 with photo- 
spheric colours, the colour temperatures range from ~ 1200K for 
NGC1978 IR1 to around 400 K for the coolest objects. 



6 CONTINUUM COLOURS AND OPTICAL DEPTH 
6.1 Manchester colours 

The IRS colour-colour diagram based on the Manchester bands for 
our targets is shown in Fig. [5] Error bars are indicated where the 
size of the bars exceeds the size of the symbol. We note that the 
errors are calculated from the tabulated noise on the spectrum, and 
do not include systematic effects such as calibration uncertainties 
or background subtraction. The majority of the stars fall on a well 
determined sequence. There three exceptions: two stars are much 
bluer in [6.4] — [9.3], and one star is too red in [16.5] — [21.5] by 
approximately 0.25 mag. 

The first two stars are GRRV38 and NGC 1978 IR4. Both 



of these have spectral energy distributions largely consistent with 
naked carbon stars. They are also faint in the IRS wavelength range 
and the colours have larger uncertainties than for the other stars, 
but the errorbars shown on the figure indicate that the flux uncer- 
tainties do not affect their positions in this diagram. The last star 
is MSX LMC 663, in which the [21.5] band includes an additional 
emission component. 

All other stars show continuum emission dominated by cir- 
cumstellar dust over all wavelenghts considered here. The good 
correlation between the two colours confirms the choice of the con- 
tinuum regions. 

6.2 Comparison models 

We calculated a series of benchmark dust models, to compare the 
colours. The models used the D USTY code, with a den sity distribu- 
tion of a radiation-driven wind jElitzur & Ivezicl200lh . p ~ p' 1 - 8 . 
Models were run with an optical depth r at 1 pm ranging from 
0.01 to 100. For the central star we assumed a black body of 
T = 2800 K. The inner radius of the shell was set at a dust tem- 
perature of 1000 K. A dust mixture of 95 per cent amorphous car- 
bon and 5 per cent SiC was used. The resultant spectra were used 
to measure the flux through the bandpasses, converted to colours 
using the zeropoints of Table [5] Such models lack the important 
molecular bands and are not meant to fit the observed spectra: this 
will be the subject of a forthcoming paper. 

The models are systematically too red in [16.5] — [21.5] by 
about 0.15 mag. The drawn line in Fig.|5|shows the models, where 
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Table 6. Photometry: fluxes and colours using the Manchester narrow continuum bands for carbon stars. For the colours in the last two columns, we adopt the 
zeropoints of Tablel5l 



target F 6 . 4 F 9 . 3 F w .s F 2 i. 5 [6.4]-[9.5] [16.5]-[21.5] 

[Jy] [Jy] [Jy] [Jy] [mag] [mag] 



TRM72 





211±0.001 





201 ±0.000 





119±0.001 





079±0.001 





758±0.005 





128±0.015 


IRAS 04557 





133±0.001 





174±0.001 





128±0.001 





102±0.001 


1 


104±0.010 





336±0.009 


IRAS 05009 





228±0.001 





247±0.001 





156±0.001 





112±0.001 





898±0.005 





220±0.013 


IRAS 05112 





231 ±0.002 





309±0.001 





233±0.001 





192±0.001 


1 


127±0.008 





336±0.007 


IRAS 05113 





160±0.001 





209±0.001 





154±0.001 





115±0.001 


1 


098±0.008 





256±0.011 


IRAS 05132 





143±0.001 





181±0.001 





132±0.001 





098±0.001 


1 


068±0.008 





247±0.016 


IRAS 05190 





225±0.002 





336±0.000 





295±0.001 





242±0.001 


1 


247±0.010 





363±0.007 


IRAS 05278 





604±0.003 





807±0.002 





563±0.004 





454±0.002 


1 


126±0.006 





341 ±0.008 


GRRV 38 





067±0.001 





037±0.001 





021 ±0.001 





014±0.000 





165±0.033 





154±0.032 


IRAS 05295 





150±0.001 





171±0.000 





115±0.001 





090±0.001 





953±0.006 





306±0.009 


IRAS 05360 





107±0.001 





135±0.001 





089±0.001 





074±0.001 


1 


062±0.006 





367±0.014 


MSXLMC 219 





112±0.001 





148±0.001 





097±0.000 





080±0.001 


1 


1 12±0.009 





370±0.009 


MSXLMC 341 





109±0.000 





126±0.001 





078±0.001 





060±0.000 





963±0.006 





291±0.012 


MSX LMC 349 





141±0.001 





166±0.000 





103±0.001 





080±0.001 





994±0.007 





300±0.009 


MSXLMC 441 





137±0.001 





181±0.001 





119±0.001 





096±0.000 


1 


115±0.013 





335±0.008 


MSX LMC 443 





083±0.001 





092±0.000 





055±0.000 





042±0.000 





918±0.008 





302±0.014 


MSX LMC 494 





094±0.000 





073±0.000 





035±0.000 





024±0.000 





540±0.006 





170±0.018 


MSXLMC 601 





115±0.001 





124±0.000 





074±0.000 





052±0.001 





894±0.006 





191±0.012 


MSX LMC 663 





116±0.001 





102±0.000 





078±0.000 





072±0.000 





677±0.006 





478±0.007 


MSX LMC 679 





135±0.001 





144±0.000 





085±0.001 





069±0.000 





886±0.006 





346±0.010 


MSX LMC 743 





093±0.000 





113±0.000 





067±0.000 





053±0.000 


1 


026±0.006 





311±0.009 


MSX LMC 749 





238±0.001 





196±0.000 





101±0.000 





068±0.001 





602±0.003 





152±0.011 


MSX LMC 754 





086±0.000 





072±0.000 





033±0.000 





023±0.000 





614±0.009 





184±0.014 


MSX LMC 967 





182±0.000 





144±0.001 





074±0.001 





051 ±0.000 





558±0.005 





180±0.011 


NGC 1978 MIR1 





112±0.001 





1 16±0.000 





072±0.001 





057±0.000 





852±0.006 





318±0.009 


NGC 1978 IR4 





042±0.000 





021±0.000 





012±0.000 





007±0.000 





079±0.024 





008±0.031 


NGC 1978 IR2 





045±0.001 





022±0.000 





012±0.000 





006±0.000 





061±0.025 


-0 


134±0.026 


NGC 1978 IR1 





106±0.000 





084±0.001 





039±0.000 





025±0.000 





554±0.007 





110±0.016 


TRM88 





098±0.001 





094±0.001 





052±0.001 





037±0.000 





769±0.009 





196±0.013 



we shifted the model tracks down by 0.15 mag to compensate for 
the offset. The cause of the offset is not known. 

The models for r = 0.01 yield a colour of [6.4]-[9.3]= 0.56. 
At these wavelengths, a 2800 K black body is considerably red- 
der than a Rayleigh- Jeans curve. The observed cutoff in this colour 
(Fig- 13 is consistent with the models, but the two stars with pho- 
tospheric spectra are considerably bluer. This blueshift reflects the 
effect of photospheric molecular bands, such as suppression of the 
[9.3] flux by the proposed C3 band. In support of this supression 
by C 3 , the [6.4]-[9.3] colour of the symbiotic C star MSX LMC 
663, which also has little dust but different photospheric molecular 
bands, is in good agreement with the models. 

We conclude that [6.4] — [9.3]< 0.5 is indicative of a naked 
carbon star. The implication is that the Manchester [9.3] band may 
not define a stellar continuum for a naked carbon star, but it can be 
used instead to judge the strength of the absorption band. 



6.3 Optical depth 

High optical depth in a circumstellar shell affects the spectral en- 
ergy distribution in two ways. Extinction within the wavelength of 
interest suppresses the bluest bands, and absorption of optical and 




0.5 1 9 8 7 6 5 

[6.4] - [9.3] [9.3] 



Figure 10. The K magnitude versus the [6.4]— [9.3] colour and versus the 
[9.3] magnitude. Symbols are as defined in Fig.BI squares and triangles in- 
dicate stars with both SiC and MgS. The drawn lines refer to the benchmark 
models of Section l6*2l 

near-infrared stellar emission lowers the dust temperature further 
out. 1 This is illustrated for our sample in Fig. llQI 

The left panel shows the K-magnitude versus the [6.4] — [9.3] 
colour. The relation between redder colour and fainter K-band mag- 
nitude in Fig. llOl confirms the effects of optical depth. The right 
panel shows the K-band versus the [9.3] band magnitude: the lat- 

1 In individual cases, density variations which differ from an r~ 2 distribu- 
tion can lead to dust temperatures which mimic optical-depth effects. 
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oxygen stars 




Figure 11. Source classification diagram based on lEgan et alj 120011) . The 
hashed area in the left panel indicat es the colour range of m odel carbon 
stars, taken from the SKY model of Wainscoat et ah 1 1992). The hashed 
region in the right panel shows the same for oxygen-rich stars. The sym- 
bols are our carbon-star sample. The dashed line in the left panel indicates 
the observed locus of obscured carbon stars. The dash-dotted line in the 
right panel shows the colour sequence for oxygen-rich stars, as derived from 
DUSTY models. 



ter is not affected by extinction within the range of optical depths 
probed by our sample. It is also close to the peak of the dust spec- 
tral energy distribution, and varies little as function of tempera- 
ture. Indeed, the distribution of [9.3] magnitudes remains fairly nar- 
row while the K-band magnitude fades. (One outlier in this plot is 
IRAS 05278—6942, which has a higher bolometric luminosity.) At 
brighter K, the [9.3] magnitude becomes a little fainter, tracing the 
effect of less dust. 

The drawn curves show the same models as discussed in sec- 
tion |^| For a single luminosity, it shows the expected behaviour 
with increasing optical depth, excluding the effects of molecular 
bands. As the depth of the shell increases, first the [9.3] magnitude 
becomes brighter, before the K-band begins to fade sharply beyond 
a critical optical depth. The curves follow the general behaviour of 
the observed spectra, but we again note that no attempt has been 
made to fit the data. 

The optical depth is a measure of the mass-loss rate. The mod- 
els imply that the K — [9.3] colour can be used as an indicator for 
the optical depth. The dashed line in the right panel of Fig llOl shows 
a constant K — [9.3] corresponding to the onset of K-band fading. 
The filled symbols (squares and triangles) show sources with 30- 
/j,m MgS emission. The dashed line shows that the presence of es- 
pecially this band is related to the optical depth in the shell. 



7 DISCUSSION 

7.1 Carbon star colour selection 

Almost all sources presented here are carbon stars. Only two 
oxygen-rich objects were found in our entire sample of LMC stars, 
and only one is classified as an AGB star. 

On one hand, this result is entirely expected. The stars at the 
beginning of the mass-loss sequence shown in Fig^a) are already 
optical carbon stars and it is expected that more evolved stars on 
the sequence would also be C stars. Such an expectation would 
only break down for massive and relatively rare AGB stars that lie 
above the the sequence of Fig. 0a) since such stars do not ever 
become carbon stars. 

On the othe r hand, based on the colour classification of 
lEgan et ail 1200 lh . a higher fraction of oxygen-rich dust shells 
would be expected. The colour criterion to separate oxygen-rich 



and carbon-rich stars, suggested by E ganetalji200lh . is based on 
the J— K versus K— A colours. The diagram is shown in Fig. II II 
The hashed regions indicate the colour range shown by their model 
carbon stars (left panel) and oxygen star s (right panel), based on 
the SKY model of lWainscoat et alJ Jl992h . Overplotted are the lo- 
cation of the stars in our sample. These mostly have colours redder 
than shown by any of the model carbon stars included in the SKY 
model. The distribution of our stars suggests that the SKY model 
can be extrapolated to redder stars. For obscured stars, there is con- 
siderable overlap in continuum colours between oxygen-rich and 
carbon-rich stars. 

There is however an effect from the silicate feature wich domi- 
nates the 10-/^m spectrum of oxygen-rich stars. The effect is shown 
for the IRAS 12-/im band bv lvan Loon et alj ll997h . At low mass- 
loss rates, the silicate is strongly in emission and enhances the 
K— [12]. At higher mass-loss rates, the silicate becomes a strong 
absorption feature which reduces the K— [12]. At a particular J— K, 
oxygen-rich stars with low mass-loss rates are redder in K— [12] 
than carbon stars, whilst for high mass-loss rates the opposite is 
true. Inevitably, there is a range where classification is not possible 
using these colours. 

The K— A colour shows a similar intrusion by silicate, but the 
MSX A-band (6.8-10.8 /im) contains more continuum and less sil- 
icate than does the IRAS 12-/im band, and therefore the flux in 
the A-band does not decline as much when silicate goes into ab- 
sorption. We calculated similar models as in Section l6~2l but for 
ox ygen-rich dust. The stel lar temperature was taken as 3200 K, and 
an Ossen kopf etal] ll992h warm 'circumstellar' silicate model was 
assumed. (These silicate optical constants may not be optimal for 
circumstellar dust.) The dashed line in the right panel of Fig. 1111 
shows the resultin g track with optical depth. It does indeed traverse 
the region used bv lEgan et alJ y001) but the K— A colour reaches 
an effective limit whilst the J— K colour continues to redden. The 
accuracy of this track is limited. The carbon-rich models of Sec- 
tion l6.2l could not be used, because the molecular bands are rather 
important for the stellar and circumstellar spectra at these wave- 
lengths, and these bands are missing from the DUSTY models. In- 
stead we draw a line through the observed data points, which are 
known to be true locations of carbon stars. At the highest mass-loss 
rates, this track is likely to steepen considerably and to run close to 
the model-track for oxygen-rich stars. At low mass-loss rate C and 
O-rich stars are well-separated, but not at high M. 

We used ISO/SWS and 2MASS data of Galactic stars to fur- 
ther examine the usefulness of n ature of the J— K vs K— A diagram. 
The ISO/SWS spectra l Kraeme r et alj 120021) are convolved with 
the MSX A-band filter transmissio n curv e. The zero magnitude for 
MSX is taken from ICohen eTai] feOOOh . The ISO/SWS stars are 
cross-correlated with the 2MASS data base and SIMBAD for iden- 
tifications. Selection criteria are 50 Jy < F(A) < 10000 Jy for 
ISO/SWS, and 2MASS flags either a,b,c,d in all of the JHK bands. 
This restricts the Galactic sample to 271 spectra. Some stars were 
observed several times and all of the A-band magnitudes were anal- 
ysed individually. Spectral types range from B-type to M- and C- 
types, and post-AGB stars and PNe (all of them are WC), as well 
as young stellar objects are included. The M-type stars are all gi- 
ants, sub-super giants, or super giants; no dwarfs or sub-dwarfs are 
selected. 

The J— K vs K— A diagram for this Galactic sample is shown 
in Fig. 1121 Early-type stars have increasing J— K while K— A re- 
mains constant, due to continuum H~ absorption. This H~ ab- 
sorption is still important for early M-type stars, while circum- 
stellar excess affects both colours among late M-type stars. C-type 
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Figure 12. The J— K vs K— A diagram for evolved stars observed by ISO 



refers to spectroscopically known carbon stars; these form a rela- 
tively narrow sequence. The isolated carbon star at [K— A]~ 15 is 
AFGL 5625: it shows a feature-less spectrum (i.e. no SiC), except 
for 14-^im C2H2. The OH/IR stars show a large range of colours; 
OH/IR stars with [K— A]~ 4—10 show silicate emission with self- 
absorption on top (e.g., CRL 2199, WX Psc), whilst OH/IR stars 
with [K— A]> 10 show silicate absorption. 

The separation between C-type and M-type stars with low 
mass-loss rates (2 < K — A < 3) is due to photospheric, 
near-infrared molecular bands. For oxygen-rich Mira stars, the 
Ks-flux is suppressed by water absorption (e.g.. iTei et al. 2003; 
iMatsuura et aljEo02bh . For carb on stars, the J-ba nd flux is sup- 
pressed by C2 and CN absorption jLoidl et all200ll) . At high mass- 
loss rates, the carbon stars have well-defined colours which over- 
lap with the (more scattered) oxygen-rich s tars. Ex tremely red stars 
have no counterpart in the LMC sample of lEgan e t al. 1 2001 ). This 
may be because of the 2MASS detection limit. 

At high mass-loss rates, the Egan et al. diagram cannot dis- 
tinguish carbon and oxygen-rich stars on an individual basis. How- 
ever, stars on the carbon-star sequence are statistically more likely 
to be carbon stars. Our selected stars fall on this sequence. At low 
mass-loss rates the Egan et al. colours of our targets are consistent 
with a carbon-rich nature. The sole oxygen-rich AGB star, IRAS 
05003—6712, has colours consistent with the carbon-rich sequence 
(J-K= 2.95, K-A= 3.8), but at the blue end of the range where 
the two types overlap. 



7.2 Molecules 

All molecules suggested to be present (C3, HCN, C2H2) are so- 
called parent molecules, which form in the photosphere, and are 
found throughout t he outflow, up to the p hotodissociation region 
at » 1-3 10 16 cm jMillar & Herbsilll994[ iMillar et alJl200ll) . Es- 
pecially C2H2 is fundamental in the chemistry leading to large or- 
ganic molecules such as PAHs. The abundance of these molecules 
is an important para meter for our un derstanding of chemistry in 
low-metallicity shells, van Loon et alj|1999a) noticed that the 3.1- 
/im band in LMC stars was of the same strength as in Galactic 
stars, and suggested that the lower metallicity was compensated by 
a higher C/O ratio. lMatsuura et aljl2002al l2005) find that C 2 H 2 is 
more abundant in LMC and SMC stars than in Galactic compari- 
son stars, based on the 3.1 and 3.%-{im absorption bands. Chemi- 
cal models 1 Mat suura et all2002al) confirm the sensitivity of C2H2 
to the C/O ratio. The carbon enhancement from third dredge-up 
leads to a higher C/O ratio for lower initial oxygen abundance. 
Ivan Loon et alj 120061) confirm that HCN is weak compared to 
C2H2 and that CS is absent, using the shape and structure in the 
3.8-/im absorption. They also show that the bands become sharper 
for stars with stronger pulsation, and interpret this as evidence for 
cold molecular gas at some distance above the stellar photosphere. 

The equivalent widths of the molecular bands are listed in Ta- 
ble|4| Fig. ll3| plots these values for the C2H2 7.5 fim and 13.7 
bands as a function of the [6.4] — [9.3] colour. The squares in the 
panels indicate the LMC stars. Especially the 13.7-/im band shows 
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Figure 13. The strength of the C2H2 features (7.5/xm and 13.7^m), as 
function of the [6.4]-[9.3] colour. Symbols are the same as in Fig.rm 



a narrow range in equivalent width. The 7.5-^im band show a large 
range, but neither correlates well with infrared colour. The two 
'naked' carbon stars show stronger bands. The equivalent width 
of the 13.7-^im feature excludes the wings of the broad component, 
and may therefore underestimate its strength for hot (photospheric) 
gas. 

Sloan et al . (2006) compare the strengths of these features for 
Galactic carbon stars and SMC carbon stars. They find that for 
stars with [6.4] — [9.3] > 0.6, the C2H2 bands are stronger in the 
SMC than in the Galaxy, the difference being more readily appar- 
ent at 7.5 fim than at 13.7 fim. The difference of strength between 
Galactic, LMC and SMC stars is interpreted as due to higher abun- 
dances of C2H2 in th e SMC, in accordance with the prediction of 
iMatsuura et alj jJoOSV Fig.[B]includes the Galactic and SMC sam- 
ples. The equivalent widths are measured in the same way for all 
stars. The Galactic sample (based on ISO spectra) is divided in 
'blue' and 'red' stars, based on whether the [6.4] — [9.3] colour is 
bluer than the lowest extinction DUSTY model (Section[62) : blue- 
ward colours are indicative of naked carbon stars. 

The Galactic sample has a large range of [6.4] — [9.3] colour. 
The SMC sample in contrast shows a very limited colour range. For 
the same mass-loss rate one may expect lower metallicity stars to 
have lower optical depths because of a lower dust abundance. The 
SMC stars should show a lower [6.4]-[9.3] than the LMC stars, 
whilst the Galactic stars should show the reddest colours. This ef- 
fect is best seen for [6.4]-[9.3]> 0.9. In the models of Section l6~2l 
an increase of a factor of 3 in optical depth reddens the [6.4] — [9.3] 
from 0.8 to 0.9, from 1.05 to 1.3, and from 1.3 to 1.8. The metallic- 
ity effect on colour is therefore much stronger for the redder stars, 
consistent with the distribution in the figure. 



The bluest objects, which are expected to be naked carbon 
stars, are almost exclusively Galactic stars, with two LMC objects. 
This is a selection effect, in that only mass-losing stars were se- 
lected for the various Magellanic Cloud surveys. For the redder 
stars, the shift in colour with optical depth suggests the Galactic 
stars with [6.4] — [9.3] between 1.2 and 1.8 should be compared 
with the LMC stars between 1.0 and 1.2, and the SMC stars around 
0.9. In the top panel, this shows a significant relation, with the 7.5- 
fim feature becoming much stronger for lower metallicity. For the 
bluer mass-losing stars, where there is little change of colour with 
optical depth, the LMC stars have a stronger 7.5-fim band than 
Galactic stars, but there is no clear offset between LMC and SMC 
stars. 

Overall, this provides some evidence for the suggestion that 
the C2H2 bands are stronger in lower metallicity environments. For 
the 13.5-/im band, a small shift is also seen between the Galactic 
and the LMC stars, although less convincing than at 7.5/im, but the 
LMC and SMC stars show comparable equivalent widths. 

The fact that the C2H2 bands show a relatively small range in 
equivalent width for a large range of colour temperatures, confirms 
that the molecule is located throughout the shell. 



7.3 Dust 

7.3.1 Temperature dependence 

The Spitzer spectra show that SiC dust bands are observed for most 
or all carbon stars, but the MgS band is present in only a subset. 
This already indicates that the MgS mineral forms under more re- 
stricted conditions. The occurrence and strength of the two bands is 
plotted for the LMC stars in Fig. 1141 as function of the continuum 
temperature (derived from [16.5] — [21 .5]) and of the [6.4]-[9.3] 
colour, as tracer of optical depth. The sources are put in three dif- 
ferent groups: the circles indicate stars without an MgS feature; tri- 
angles show stars with MgS and with weak SiC (defined as a line- 
to-continuum ratio less than 0.14); squares show stars with MgS 
and stronger SiC. 

The figure shows that the presence of MgS correlates strongly 
with the dust temperature. There is a sharp division at a (black 
body) temperature of T = 650 K: none of the stars with hot- 
ter dust show MgS, whilst all stars with cooler dust do. Between 
dust temperatures of 500 and 600 K the MgS feature is relatively 
weak, and at lower temperatures the feature is stronger and shows 
a large range of line-to-continuum ratio. In contrast, SiC tends to 
be weaker at high dust temperatures but is almost always present. 
Panel d indicates that increasing SiC strength is accompanied by 
rapidly increasing MgS. The difference can be understood in terms 
of the formation mechanism of the different dust species. Solid SiC 
condenses directly from the gas, at high temperatures. MgS grows 
as a surface on pre-existing grains iNuth et alll985l) : this process 
only starts at temperatures around 600 K, and is complete around 
300 K. 

Fig. 1 151 shows a less common way to look at MgS, using the 
relation with pulsation period of the star (taken from Table|2j. The 
double-period star is not plotted. The MgS feature is seen in gen- 
eral for stars with periods longer than 650 days. The relation is an 
indirect one, in the sense that the longest period stars tend to have 
the highest mass-loss rates. 
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Figure 14. The strength of the 11-micron SiC and the 30-micron MgS bands, as function of the dust continuum temperature (derived from [16.5]— [21.5]), 
and the [6.4]— [9.3] colour. Circles indicate stars without an MgS feature; Squares: stars with an MgS feature and a strong SiC feature; Triangles: stars with 
an MgS feature and a weak SiC feature. 



7.3.2 Metallicity dependence 

Fig. ll6l compares the strengths of the SiC and MgS fe atures in our 
LMC sample with the SMC and Galactic samples of Slo an et alj 
(2006). Sloan et al. find that the SiC and MgS features are weaker 
in the SMC than in the Galaxy. This is expected, as Mg, Si and S 
are not produced in AGB stars and their abundance should scale di- 
rectly with metallicity. Our data show a more complicated picture, 
however. 

The line-to-continuum ratios are plotted against [6.4] — [9.3], 
which is a measure of optical depth. As before, for the higher-mass- 
loss stars the optical depth is less in the SMC, while at low optical 
depth SMC, LMC and Galactic stars can be compared at the same 
value for [6.4] — [9.3]. The top panel shows that at low optical depth 
([6.4] — [9.3] ~ 0.5), there is a clear sequence with the SiC line-to- 
continuum ratio diminishing with decreasing metallicity. At higher 



optical depth, the three samples indicate very similar ratios. Stars 
with 'naked carbon star' colours show a high line-to-continuum 
ratio, related to the weakness of the continuum, with a sequence 
clearly offset from the dusty stars. It is possible that here the con- 
tinuum used to measure the band strength is affected by molecular 
absorption bands, giving an artifically enhanced feature. The SiC 
band is flanked by two photospheric bands, attributed to C2H2 and 
C3, and a true continuum cannot easily be defined. 

The MgS feature is illustrated in the bottom panel of Fig. 1 161 
There is an indication that the MgS band is weaker in the SMC than 
in comparable Galactic and LMC stars. The LMC stars show a very 
large scatter and it is not possible to draw conclusions regarding its 
typical strength. We note that the continuum underneath the feature 
is not well defined in Spitzer spectra because of its long-wavelength 
cutoff. 

In both cases the line-to-continuum ratio in dusty stars com- 
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Figure 15. The MgS line-to-continuum ratio as function of pulsation period 
of the star. Open circles: no MgS, triangles: weak SiC; squares: strong SiC 
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Figure 16. The strength of the SiC and MgS features as a function of the 
[6.4]-[9.3] colour. 



pares the band strength with the underlying dust continuum, as- 
sumed to be due to amorphous carbon. As the carbon is dredged-up 
in the star while the mineral constituents are not, the ratio is ex- 
pected to decrease with decreasing metallicity. However, this is also 
affected by the dust condensation process. If seeds are needed to 
grow the amorphous carbon grains (e.g. TiC), the relation between 
mineral and amorphous carbon dust could be more complicated. 
Especially the SMC data indicate that the emissivity ratio between 
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Figure 17. Apparent central wavelength of the SiC feature as a function of 
the [6.4]-[9.3] colour. The dashed lin e represent s the distributions of Galac- 
tic carbon stars, according to lSloan et alj i2006). Symbols are as in 



minerals and amorphous carbon dust is lower at low metallicity, but 
the LMC data is inconclusive on this point. 

The abundance of SiC is limited by the enriched carbon and 
the unenriched silicon. The abundance of MgS is limited by two 
unenriched elements. If the condensation is 100 per cent efficient, 
both will be affected similarly by metallicity. If the condensation 
operates at lower efficiency, the formation becomes limited by the 
collision rate and the the abundance of MgS would more strongly 
dependent on metallicity than SiC ([MgSJoc Z 2 , while [SiC]tx Z, 
where Z is the fractional abundance of metals). (Temperature ef- 
fects may also contribute, as lower opacity gives rise to higher dust 
temperatures.) The spectra presented here do not support such a 
strong dependence on metallicity, but favour a linear relation. This 
supports a high condensation efficiency for MgS. Confirmation will 
require detailed modelling of the dust spectra. 



7.3.3 Central wavelengths 

The central wavelengths of the dust bands are important diagnostics 
ISpeck et all2005l:lHonv et alj|2002h of the dust environment. We 
define the central wavelength (hereinafter Asic) as the wavelength 
which disects the flux distribution in the band: i.e. after continuum 
subtraction, half the band flux arises redward and half blueward. 
This definition is more stable against noise on the spectral energy 
distribution than the wavelength of peak emission strength. How- 
ever, for MgS only part of the band is covered by the Spitzer spectra 
and the derived central wavelength reflects this. 

Fig. 1171 shows the observ ed variation of Asic, as a func- 
tion of the [6.4]-[9.3] colour. ISloan et alj 120061 have shown that 
for Galactic carbon stars with [6.4] — [9.3] >0.5, Asic increases 
monotonically with [6.4] — [9.3] colour, from ~ 11.20 /im up to 
~ 11.70 /im for the reddest source of their sample. For bluer 
colour s Asic redd ens again: their relation is indicated by the dashed 
line. ISpeck et"al] 120051) argue that red carbon-rich stars can show 
the SiC in self-absorption and full-absorption. This absorption oc- 
curs at 10.8 ^m, so that an increase of the amount of SiC dust lead- 
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ing to self-absorption will cause a redward shift of Asic- The pres- 
ence of an absorption component in all stars remains to be proven. 

The inverse shift to longer wavelength of Asic, for stars with 
[6.4] — [9.3] <0.5, is not yet explained. This colour corresponds to 
naked carbon stars, where the [9.3] band may be suppressed by a 
photospheric molecular band. This band borders the SiC feature, 
and may overlap with it in wavelength. This would have a sim - 
ilar effect as the self-absorption proposed by ISpeck et"aH i2005l) . 
and could self consistently explain why a bluer [6.4] — [9.3] colour 
(deeper 10-/im absorption) correlates with a redder SiC feature. 
Note that the broad molecular absorptions at 10 and 13 /im may 
give the impression of an intermediate emission band: the actual 
presence of SiC in such stars may be in doubt. 

In our LMC sample, we observe that Asic is redder for the 
stars with weaker SiC features (triangles in Fig. 1171 . These stars 
have slightly redder [6.4] — [9.3] as shown in Fig. 1141 This can in 
principle be interpreted as being due to self-absorption weakening 
the band, although other explanations are possible. In our sample, 
the dust continuum temperature has no discernable influence on 
Asic- 

The LMC stars with MgS (triangles and squares: cool dust) 
follow the slope of the relation of ISloan et alJ f2006), but offset 
towards redder [6.4] — [9.3]. An alternative view would be that at the 
same optical depth, LMC stars have a bluer Asic - This may reflect 
the broad C2H2 R-band at 12-14 /im, depressing the apparent red 
continuum. A stronger C2H2 band in LMC stars would shift the 
observed Asic to shorter wavelength. 

Fig. ll8l shows the central wavelength of the MgS feature (here- 
inafter AMgS, defined as for the SiC feature central wavelength) 
as a function of the temperature derived from the [1 6.5] — [21.5] 
colour . From a large sample of Galactic evolved stars, Ho nv et alJ 
1 2002) show that the central wavelength of the MgS feature is inde- 
pendent of continuum temperatures. They explain this by the fact 
that MgS is partially heated by mid-IR radiation. In our sample, 
the MgS band is observed only in stars spread over a small range 
of continuum temperatures, making it difficult to find any correla- 
tion between AMgS and the temperature. The suggestion of longer 
AivigS for higher dust temperatures is not conclusive. Our LMC 
sample shows AM g s similar to Galactic stars observed by Hony et 
al. within the same temperature range. Fig. 1191 show the shape of 
the continuum-substracted MgS feature. This figure confirms that 
there is no apparent correlation between the MgS feature shape and 
the dust continuum temperature. The short-wavelength band edge 
has a similar shape in all spectra. There is some separation between 
stars with relatively flat MgS band and stars whhere the band dimin- 
ishes towards the red, but this separation does not show an evident 
correlation with dust temperature. 



7.3.4 Oxygen-rich dust 

The current sample contains almost exclusively carbon-rich stars. 
Oxygen-rich mass-l osing stars do exist in the LM C, as found in 
the s ample of, e.g.. iTrams et al. 1 1999) (see also Diikst ra et alJ 
l2005j). Comparing th e J— K versus K— A diagram for MSX sources 
in ligan et alJ iiooil) . with the diagrams in Section l7~T1 shows that 
the majority of cocoon LMC stars have colours consistent with the 
carbon star sequence, but some stars are bluer in K— A and these 
could be oxygen-rich. 

The higher efficiency of third dredge-up at low metallicity 
turns LIMS into carbon stars relatively early on the AGB: the exis- 
tence of optical (i.e. unobscured) carbon stars show that in the LMC 
this typically occurs before the onset of the high mass-loss phase. 



30.5 



30.0 



29.5 



29.0 
300 



A - _ 



400 500 600 

Temperature [K] 



700 



Figure 18. Apparent central wavelength of the MgS feature as a function of 
the temperature (derived from the [16.5] — [21.5] colour). Symbols are as in 
Fig-IHI 



Oxygen-rich mass-losing stars will most likely be either stars with 
initial masses large than about 4M©, where hot bottom burning 
prevents the star from becoming carbon-rich, or low-mass stars. 

Dust formation in oxygen-rich stars is dependent on 
metallicity-dependent minerals, such as corundum and various sil- 
icates. Oxygen-rich stars may therefore show relatively poor dust 
formation efficiency at low metallicity. For some stars, this could 
in principle delay the superwind phase sufficiently to allow them to 
become a carbon star. Carbon stars are less affected by metallicity, 
as dust is produced from self-enriched carbon. It would be of inter- 
est to test whether the LMC stars show any evidence for a delayed 
or weaker superwind for oxygen-rich stars. 

The predominance of mass-losing carbon stars implies that 
dust enrichment of the ISM at low metallicity by AGB stars will 
be strongly dominated by amorphous carbon dust, as compared to 
the Galaxy. 

7.4 Three groups 

In Sect. I7.3.T1 we have shown that the plot of the strength of the 
dust bands versus colours or temperatures permits us to identify 
three distincts group of carbon stars in our LMC sample. The three 
groups can be defined as: (1) stars with a strong SiC feature and 
MgS, (2) stars with a weak SiC feature and MgS, and (3) stars 
without an MgS feature. 

Fig. 1201 shows the spectral shapes of stars within those three 
groups, obtained by averaging the different spectra within each 
group. 

The shape of the SiC band is essentially identical between the 
first two groups. The continuum on either side is different. The nar- 
row 13.7-/im band weakens as SiC becomes weaker, but the broad 
band between 12 and 16 ^m becomes much more prominent. This 
is best seen by comparing the slopes of the spectra. The blue slopes 
of the continua (5-9 /im) also show a difference between the three 
groups. 

Such a classifiation can be compared with that proposed by 
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iKraemer et alj l2002h . who divide dusty carbon-rich stars into red 
(CR) and blue (CE). In our sample, about 10 stars would be classi- 
fied as CE and the remainder (apa rt from the symbiotic C star) as 
CR. In the Galactic ISO sample of Kraemer et alj J2002I) . CE stars 
are somewhat more dominant. This difference is likely to be caused 
by selection effects. 



8 CONCLUSIONS 

We have presented a Spitzer survey of 29 mass-losing Asymptotic 
Giant Branch stars in the Large Magellanic Cloud. Most are classi- 
fied as carbon stars, based on the infrared spectra. One of the C stars 
is also a symbiotic star. Two oxygen-rich stars were also present in 
the original survey: one was found to be a high-mass young stellar 
object and is presented in lvan Loon et all f2005bl) . 

The spectra cover the wavelength range 5-38 /im. The short 
wavelength range, up to 16 pirn, is dominated by absorption bands 
caused by circumstellar and/or photospheric molecules. At the 
short edge, bands of CO and C3 cause a sharp drop of the stellar 
energy distribution towards the blue end of the spectral range. Ab- 
sorption bands at 7.5 and 13.7 /im are due to acetylene, C2H2. The 
13.7/im band contains both a narrow and a broad component, the 



latter tracing hotter gas. These bands are also seen in two naked 
carbon stars in our sample. 

Dust emission bands are present in our sample at 11 and at 
30 /im, due respectively to SiC and MgS. The first band is present 
in almost every object, regardless of dust temperature. MgS is only 
present when the dust temperature is below 650 K, approximately 
in half the sample. This is consistent with how these minerals form: 
SiC condenses directly out of the gas phase, but MgS is deposited 
on existing dust grains below a critical temperature. 

Because of the many bands in the spectrum, continuum def- 
inition is difficult. We define a set of four narrow bands, called 
the Manchester system, which can be used to trace the continuum 
for dusty carbon stars. The bands are centred at 6.4, 9.3, 16.5 and 
21.5 Ltm. The [6.4] — [9.3] colour is shown to be a measure of the 
optical depth in the shell. The [16.5] — [21.5] colour can be used to 
measure the dust temperature. Two naked carbon stars show much 
bluer colours in [6.4] — [9.3] than expected from models. This is 
attributed to suppression of the [9.3] band by a molecular band, 
which is apparent in many of the stars but is particularly clear 
in stars with photospheric infrared emission. A 10-/tm absorption 
band in carbon stars has previously been attributed to interstellar 
silicate, but this assignment can be ruled our for the LMC stars. A 
recent assignment to silicon-nitrite particles can also be ruled out 
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Figure 20. Spectra of the three object classes observed, the first class (top) containing objects displaying a weak SiC feature and a MgS feature, the second 
(middle) a strong SiC and a MgS feature, the last one (bottom) displaying no MgS feature. These three spectra are obtained by averaging the individual spectra 
for each class. 

The 13.7-/im band does not clearly show this. The SiC band shows 
evidence for a much lower line-to-continuum ration at the lowest 
metalicity (the SMC) but no strong difference between the LMC 
and the Galaxy. For MgS there is an indication for a sequence of 
decreasing line-to-continuum ratio with decreasing metallicity. We 
show that the central wavelength of the SiC feature dep ends on the 
[6.4] — [9.3] colour, consistent with Sl oan et all ^2006). It is how- 
ever important to realize the uncertainty of the continuum determi- 
nation, affected by the molecular absorption bands on either side. 

The increasing strength of acetylene bands with decreasing 
metallicity can be explained by the effect of carbon dredge-up: a 
fixed amount of added carbon gives a higher C/O ratio at lower 
metallicity (less O), leaving more free carbon after the formation 
of CO. This may also be the cause of the C3 bands which we see. 



for emission at photospheric temperatures. This supports the id en- 
tification with a C3 band, as proposed by Jorge nsen et alj f2000). 

The classification scheme of mass-losing stars by Ega n et all 

1200 ll) is investigated. Based on this scheme, some stars in our sam- 
ple were expected to be oxygen-rich. We show from current data, 
from dust models and from ISO spectra of Galactic stars, that for 
high mass-loss stars, the J— K, K— A colours cannot easily distin- 
guish oxygen-rich from carbon-rich stars. The fact that almost all 
our sample was found to be carbon rich reflects the fact that they lie 
on a sequence of increasing mass-loss rate, with the stars becoming 
carbon stars even before they develop thin dust shells. 

We investigate the strengths of the molecular and dust bands 
as function of metallicity, based on a comparison with carbon stars 
in the Galaxy and the SMC. The 7.5-fim C2H2 band shows evi- 
dence for increasing equivalent width with decreasing metallicity. 
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SiC and MgS do not benefit form this effect, and their abundance 
is expected to scale more directly with metallicity. 

We finally note that dust formation in oxygen-rich stars is de- 
pendent on metallicity-dependent minerals, such as corundum and 
various silicates. Oxygen-rich stars may therefore show relatively 
poor dust formation efficiency at low metallicity. Carbon stars, on 
the other hand, depend for dust on self-produced carbon and can 
have high efficiency of forming amorphous carbon dust. The ear- 
lier formation of carbon stars due to third dredge-up, together with 
a possible difference in dust formation efficiency, means that dust 
enrichment of the ISM at low metallicity by AGB stars will be dom- 
inated by amorphous carbon dust, as compared to the Galaxy. 
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